Photodisintegration of light nuclei 
for testing a correlated realistic interaction in the continuum 
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An exact calculation of the photodisintegration cross section of 3 H, 3 He and 4 He is performed using 
as interaction the correlated Argonne V18 potential, constructed within the Unitary Correlation 
Operator Method (Vucom)- Calculations are carried out using the Lorentz Integral Transform 
method in conjunction with an hyperspherical harmonics basis expansion. A comparison with other 
realistic potentials and with available experimental data is discussed. The Vucom potential leads 
to a very similar description of the cross section as the Argonne V18 interaction with the inclusion 
of the Urbana IX three-body force for photon energies 45 < (J < 120 MeV, while larger differences 
are found close to threshold. 

PACS numbers: 21.45.+V, 21.30.-x, 25.20.Dc, 27.10. +h 



I. INTRODUCTION 

One of the main challenges in theoretical nuclear 
physics is a microscopic description of the properties of 
finite nuclei, using realistic nuclear forces. For a given 
interaction model, appropriate observables showing im- 
pact on different features of the nuclear potential need 
to be tested in a many-body system through direct com- 
parison with experimental data. In few-particle systems, 
where the quantum many-body problem of nucleons can 
be solved exactly both for bound and scattering states, 
one has an optimal setting to probe off-shell properties 
of different potentials. 

For photodisintegration of light nuclei, exact calcula- 
tions predict a slightly different behavior when using dif- 
ferent two-body realistic interactions, while a more evi- 
dent effect of three-body forces is found. However, the ex- 
perimental situation is still unsatisfactory, as the data are 
not sufficiently precise to discriminate clearly among dif- 
ferent interactions models. Among light nuclei, the pho- 
toexcitation of the alpha particle has recently attracted 
much attention both in theory, where a calculation with 
realistic two- and three-body forces was carried out [l|, 
and in experiments 0, H[ i where the aim was to further 
clarify whether He exhibits a pronounced giant dipole 
resonance or not, a question that was raised in the first 
microscopic calculation of the reaction Q , where semire- 
alistic two-body forces where used. Unfortunately, exper- 
imental data do not yet lead to a unique picture. Nev- 
ertheless, a comparison of the impact on observables of 
different potentials on a theoretical level is already in- 
structive. 

These exact calculations of the photodisintegration 
cross section can also offer valuable theoretical guidance 
in a class of astrophysical studies. Nuclear absorption 
of high-energy 7-rays through excitations of giant reso- 
nances may potentially become an important diagnostic 



tool for the dense environments of compact astrophysical 
sources. For example, black holes emit an intense high- 
energy 7-ray continuum, against which nuclear absorp- 
tion features would uniquely probe the baryonic-matter 
density, and possibly some compositional information, 
near the source. Since the dependency on the potential 
model used in our calculations turns out to be small in 
comparison to the discrepancy in experiments, our pre- 
cise microscopic reaction calculations can help to inves- 
tigate interstellar gas, mostly isotopes of hydrogen and 
helium, surrounding the astrophysical 7-ray source Q. 

The purpose of this paper is to investigate the to- 
tal photodisintegration cross section of light nuclei us- 
ing the potential constructed within the Unitary Corre- 
lation Method (UCOM), denoted with Vucom- This is 
aimed at testing for the first time this interaction in a 
continuum reaction, where many disintegration channels 
are open, and simultaneously studying the capability of 
describing electromagnetic reactions via non-local two- 
body interactions versus local two-and three-body mod- 
els. Exact calculations are performed making use of the 
Lorentz Integral Transform (LIT) approach [6(, which en- 
ables to fully take into account the final state interaction, 
while circumventing the difficulties arising from the con- 
tinuum many-body scattering states. The problem is in 
fact reduced to the solution of a bound-state-like equa- 
tion, which we solve with an hyperspherical harmonics 
(HH) basis. 

The paper is organized as follows. In Sec. [TT| the the- 
oretical background is set: an overview of the UCOM 
is given, then the LIT method is briefly summarized. 
In Sec. IIIII results are shown and conclusions are finally 
drawn in Sec. IIVI 



II. THEORETICAL OVERVIEW 



A. The Unitary Correlation Operator Method 
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Different models for the nucleon-nucleon (NN) interac- 
tion are found in the literature, as e.g. the Argonne V18 
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potential (AV18) 0, the CD Bonn [| and the Nijmegen 
[9f potentials. These modern two-body interactions re- 
produce the experimental NN data with high precision, 
but they undcrbind nuclei with A > 3. This drawback 
has traditionally been overcome introducing phenomeno- 
logical three-body forces fitted to reproduce experimen- 
tal binding energies and first excited states of light nu- 
clei. Furthermore, chiral perturbation theory can provide 
a systematic way to build two-, three- and more-body 
forces [13, EH- The above mentioned potential models, 
including two- and three-body forces have been used to 
calculate bound state properties of light nuclei (up to 
A = 12 and more), mainly within the Green's function 
Monte Carlo (GFMC) method [llJllEland the no- 
core shell model (NCSM) approach lHT 1 6L 1 1 7j] - However, 
the addition of a three nucleon potential requires a very 
intensive computational effort when used in many-body 
systems. In Ref. [l8| it was shown that different, but 
phase equivalent two-body interactions are related by a 
unitary non-local transformation. It is expected that the 
addition of non-locality to the NN interaction could re- 
duce or even cancel the need of three-body force, leading 
to a big simplification of many-body calculations. Based 
on this idea, a new category of NN potentials has then 
emerged, among them one can recall: (i) the Inside Non- 
local Outside Yukawa (INOY) interaction by Doleshall 
et al. 1 1 9t | , (ii) the J-matrix Inverse Scattering Potential 
(JISP) by Shirokov et al. [H SH and (iii) the Unitary 
Correlation Operator Method (UCOM) by Feldmeier et 
al. [H, S3, Q. The Vucom potential derived fr om the 
AV18 interaction has already been used as universal in- 
put in quite a variety of many-body techniques, ranging 
form Hartree-Fock calculations 25], to Random Phase 
Approximation [26| and to nuclear structure calculations 
in the framework of Fcrmionic Molecular Dynamics, see 

e.g. nam. 

In the UCOM framework the dominant short range- 
correlations, induced by the repulsive core and by the 
tensor part of the NN interaction, are explicitly described 
by a state independent unitary transformation. When 
the correlation operator is applied to the Hamiltonian, a 
phase-shift equivalent correlated interaction is obtained. 
A correlated operator is defined via a similarity transfor- 
mation 



body contributions, n — 1, . . . , A, 
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where 0^ n > indicates the irreducible n-body part. In the 
UCOM one usually employs a two-body cluster approx- 
imation, neglecting CH 3 ! and higher cluster orders. The 
effect of this approximation was tested in the framework 
of NCSM, where one can perform exact few-body calcu- 
lations 29]. In particular, in this reference it was shown 
that the omitted three- and more-body terms of the clus- 
ter expansion can be tuned changing the range of the ten- 
sor correlators in order to compensate the genuine three- 
body force to a large extent, still preserving the phase 
shift equivalence. The value of the tensor "correlation 
volume" 4 1,0) = 0.09 fm 3 (for details see ^) was found 
to give the best description of binding energies of 3 H and 
He on the Tjon line. In this paper we will use this very 
same value of Ig 1 ' ^ in order to investigate whether the 
non-local two-body Vucom interaction, that minimizes 
the effect of three-body forces on the binding energies, 
can also describe the continuum photoabsorption cross 
section. 

As already pointed out in [l8j , phenomenological non- 
local interaction terms introduce modifications on the 
electromagnetic current operator, conserved by gauge in- 
variance. In fact, if the nuclear potential V does not 
commute with the charge operator, then two-body cur- 
rents, usually called Meson Exchange Currents (MEC), 
have to be introduced, for consistency. In case of phe- 
nomenological non-local potentials the explicit construc- 
tion of a consistent MEC could be rather involved. The 
study of electromagnetic reactions at low energies, where 
one can use the Siegert theorem (see e.g. [30(), allows to 
investigate the role of implicit electromagnetic exchange 
mechanisms without requiring an explicit knowledge of 
the two-body current operator. For this reason, the test 
of a given potential model on the prediction for photoab- 
sorption cross section is very important. 



B. The Lorentz Integral Transform approach 

The total photoabsorption cross section is given by 

a(iv) = 4n 2 aujR(uj) , (3) 



where C r and Cq are the unitary radial and tensor cor- 
relation operators, respectively. The correlator C r intro- 
duces a radial distance-dependent shift which keeps nu- 
cleons away from each other when their uncorrelated dis- 
tance is smaller than the range of the strongly repulsive 
core of the NN interaction. Cq produces tensor correla- 
tions by further spatial shifts perpendicular to the radial 
direction, depending on the orientation of the spins of 
the two nucleons with respect to their distance. 

In a many-body system the correlated operator of 
Eq. {1} becomes an A-body operator with irreducible n- 



where a is the electromagnetic coupling constant and, at 
low photon energy uj, -R(w) is the inclusive unpolarizcd 
dipole response function, generally defined as: 



Here, Jo and Mq indicate the total angular momentum of 
the nucleus in its initial ground state and its projection, 
while | Wo//) and E Q /f denote wave function and energies 
of the ground and final states, respectively. The dipole 
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operator is 



(5) 



where Zi and rf are the third components of the position 
in the center of mass reference frame and isospin of the 
z-th particle, respectively. The dipole approximation has 
been proven to be very good at low photon energy for 
the deuteron [3l[ and for the triton case 32]. With the 
dipole operator the major part of the meson exchange 
currents (MEC) is implicitly taken into account, via the 
Siegert theorem. 

In the LIT method @ one obtains R(u>) after the in- 
version of an integral transform with a Lorentzian kernel 



L((TR,cri) = du> 



(uj - cr R ) 2 + aj 



WW) 



(6) 



where W) is the unique solution of the inhomogeneous 
"Schrodinger-like" equation 

0-Eo-<rn + uri)\*) = D x \*o). (7) 

Because of the presence of an imaginary part 07 in Eq. (|7|) 
and the fact that the right-hand side of this equation is 
a localized state, one has an asymptotic boundary con- 
dition similar to a bound state. Thus, one can apply 
bound-state techniques for its solution. The response 
function R(lv) is recovered via an inversion of the inte- 
gral transform of Eq.© (for inversion methods see Refs. 
[33| and [Uj]). Up to now this method has allowed cal- 
culations of electromagnetic reactions for 3 H and 3 He 
[32l . l35l I3H [37] ] and recently even for 4 He [l| with realis- 
tic two- and three-body forces. Nuclei with six- [H, [39[ 
and even seven nucleons [4(| have been investigated 
with semirealistic interactions. A similar formalism has 
been applied to describe exclusive electromagnetic pro- 
cesses in the four-body system with semirealistic poten- 
tials [H, [43|]. Recently, the method has been also 
applied to the harmonic oscillator basis in the framework 
of NCSM Q. 

In the present paper we calculate the LIT using an 
HH expansion of the internal wave function. For the 
antisymmetrization of the wave function we make use of 
the powerful algorithm developed by Barnea et al. |45l 
HE E3 ■ The HH approach, which is natural in coordinate 
space for a local interaction, has been recently extended 
to the use of non-local potentials given in terms of two- 
body matrix elements between harmonic oscillator (HO) 
states (48|, like the JISP potential. In this paper we 
will make use of the same method using as input the 
correlated HO matrix elements of the AV18 potential, 
constructed within the UCOM. 
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FIG. 1: (Color online) Ground-state energies of 4 He and 3 H 
with the Vucom potential for different values of the HO pa- 
rameter hfl as a function of the HO excitations N max allowed 
in the NCSM model space and of the maximal value of the 
HH grand-angular momentum quantum number K max in the 
HH expansion (n max = 15). 



AV18 potential, the correlated Vucom contains deriva- 
tive operators. In principle, one could extend the HH for- 
malism to the use of such terms and consider the rather 
involved operator form of the Vucom, the result being 
a more intensive numerical computation for each matrix 
element. The alternative approach we use, consists in 
representing the two-body potential operator on an HO 
basis, similarly to the case of the JISP potential. The 
input of our HH calculation are the correlated relative 
two-body HO matrix elements, 
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where n,£, s and j are radial quantum number, orbital 
angular momentum, spin and total angular momentum 
of the two-body sub-system, respectively (isospin is omit- 
ted for the sake of simplicity). In the operator form the 
UCOM is independent on the HO parameters, therefore 
the HO representation should be viewed as a parameteri- 
zation rather than a formulation of the potential. In fact, 
we expand the two-body potential as 



III. RESULTS 

We begin the discussion with the results for the ground 
state energies of 4 He and 3 H calculated in the HH expan- 
sion with Vucom ■ Even if the starting point is the local 
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The sum over two-body quantum numbers n, n' and £, £' 
has to be performed up to maximal value of n max and 
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TABLE I: Ground state energies in MeV with the Vucom po- 
tential. NCSM results from [2S|. HH calculations performed 
with n ma x=15- 



Method 


4 He 


3 H 


NCSM 


-28.4(1) 


-8.32(3) 


HH 


-28.57(3) 


-8.27(2) 


Nature 


-28.30 


-8.48 



(■max for the radial quantum number and angular mo- 
mentum, respectively. The expansion has to be pushed 
forward, till independence of observables, like binding en- 
ergies, on the HO parameters is reached. We find out 
that good convergence of the angular part is reached with 
(■max = 6, while at least n max = 15 oscillator quanta are 
needed for the radial part. 

Since we work in an hyperspherical harmonics many- 
body Hilbert space, convergence as a function of the 
hyperspherical grand angular momentum K max needs 
to be further investigated. In Fig. [TJ we compare the 
HH results versus K max with the NCSM data from [2^| 
as a function of the HO excitation allowed N max = 
J2i 2 n i + (%i where i runs over the A — 1 Jacobi coor- 
dinates. Results for different values of the HO param- 
eter Ml are presented. As one can note, the NCSM 
^-dependence is rather strong, while the HH approach 
is almost HO parameter free. The small residual /in- 
dependence, of the order of 0.7%, is due to the trun- 
cation in Eq. ©. The HH convergence as a function 
of K max is much faster than the NCSM convergence in 
N m ax- This is strongly manifested in case of the more 
extended nucleus of 3 H, where convergence of the NCSM 
is reached only with N max — 48. The reason for that lies 
in the fact that in the HH-approach matrix elements with 
oscillator quanta up to n max = 15 of the two-body sub- 
system are considered for each K max value, and not only 
beyond N max = 30, as for the NCSM. In Tabled the re- 
sult for ground state energies of Vucom for 4 He and 3 H 
are summarized. The NCSM results correspond to the 
Ml that yields the minimal energy for the largest model 
space used, whereas for the HH results the mean value of 
the different Ml has been taken. The two methods agree 
with each other. The slightly lower 4 He and the slightly 
higher 3 H binding energies in HH are related to the dif- 
ferent Hilbert spaces. They are consistent with the fact 
that in the HH-approach with ( max = 6 and n max = 15 
one takes into account higher two-body HO excitations 
than with N max = 16, but less than with N max = 48. 

The ground state wave function \^o) and energy Eq 
are then used in Eq. ([7]) to address the problem of the 
continuum with the LIT method. In Fig. [21 we show the 
transforms L(an,ai — 20 MeV) in case of 4 He as a func- 
tion of the parameter for Ml fixed to 20 MeV. The 
convergence in terms of the HH expansion is studied as 
a function of the K max : the even/odd value of K max is 
due to parity difference of the expanded states, j^o) and 
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FIG. 2: (Color online) 4 He: (upper panel) convergence of the 
LIT, with fixed 07 = 20 MeV and HO parameter Ml = 20 
MeV, for different values of the hyperspherical grand-angular 
momentum quantum number K max in the HH expansion 
(n m ax = 15); (lower panel) convergence of the photoabsorp- 
tion cross section. 



D^o), respectively. After inverting the various trans- 
form and making use of Eq. ([3]) one can observe the con- 
vergence of the photoabsorption cross section. Full con- 
vergence is reached with K max = 14/15 both for the LIT 
and for a~ r Similar behavior in K maxi though slightly 
weaker, is obtained for the Ml values of 40 and 60 MeV. 
In an analogous way we have investigated the 3 H case, 
where full convergence is reached with K max = 16/17. 
The broader structure of 3 H with respect to 4 He makes 
the rate of convergence slightly slower. 

At this point it is convenient to compare the converged 
cross sections obtained with different HO parameters to 
investigate the Ml dependence in the continuum. As one 
can see in Fig. [3] in case of 4 He one gets stable results 
for (7 7 with n max = 15. The residual Ml dependence can 
be interpreted as the numerical error of our calculations: 
the indetermination is maximally 1.5% in the dipole res- 
onance peak region and never exceeds 6% for higher en- 
ergies. On the contrary, in case of 3 H the HO parameter 
dependence is still rather strong with n max = 15, thus 
higher values of n max need to be considered in the ex- 
pansion of Eq. ©. In particular, we observed that for 
the larger Ml — 40 and 60 MeV values, higher n max are 
needed to obtain stable result for the cross section. To 
check the degree of convergence in the continuum one can 
investigate the first sum rules of the cross section. For 
example, we have considered the inverse energy-weighted 
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FIG. 3: (Color online) Photoabsorption cross section of 4 He 
(upper panel) and 3 H (middle panel) for different values of 
the HO parameter hfl. For 4 He a maximal value n max = 15 
in the correlated HO matrix element is used in input, while 
for 3 H n max — 25 is considered. (Lower panel) convergence 
of the sum rule m_i for 3 H as a function of n max . 



sum rule of the cross section m_i (proportional to the 
total dipole strength), evaluated as expectation value on 
the ground state [491 ] . 



V 7 (w) 



du = 



Z(Z-l) 



• 2 ) 

VP/ 

( 10 ) 

Here (r^) and (r"z p ) are the mean square point-proton 
and mean square proton-proton radii. In Fig. [3l the 
behavior of the m_i for 3 H as a function of n max (for 
converged expansion in K max ) is shown for different HO 
frequencies. One can see that for n max — 15 the discrep- 
ancy is of about 5% (1 — 2% in case of 4 He), while only 
with n max = 25 a satisfactory ^-independence (0.5%) is 
achieved. This is again due to the fact that the 3 H wave 
function has a longer range structure than 4 He, thus con- 
vergence is slower. Furthermore, the higher sensitivity of 
77i_i in n max with respect to Eq ( 3 H binding energies 
decreases of only 30 KeV going from n max = 15 to 25) 
is related to the long range nature of the operator. In 
the peak region, the final 3 H photoabsorption cross sec- 
tion presents a 1% difference between the hfl — 20 MeV 
and the hfl = 40 MeV results, and about 4% between 
the hfl = 20 MeV and the hfl = 60 MeV. The reason of 
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FIG. 4: (Color online) Photodisintegration cross section of 
triton (upper panel) and 4 He (lower panel) with different po- 
tential models: (thin line) Vucom, (dashed line) bare AV18 
and (thick line) AV18+UIX force. 



the poorer convergence for a high HO parameter in the 
continuum was already pointed out in [44| . where the 
LIT method was applied for the first time to the NCSM 
basis. Here, though we use an HH basis, the parameteri- 
zation of the potential via HO two-body matrix element 
leaves some fingerprints, i.e. a slight residual hfl depen- 
dence. As the HO potential well becomes steeper and 
steeper, the two-body wave functions in a fixed model 
space (fixed n max ) go faster to zero, consequence being 
a poorer representation of long-range operators, like the 
dipole. Thus, higher n max need to be considered to get 
the same degree of of accuracy. On the contrary, for 
smaller HO frequencies, the two-body HO eigenstates are 
closer to each other, resulting in a better sampling of the 
complex-energy continuum of the LIT, and thus in a bet- 
ter convergence of the cross section. 

In Fig. we compare the photoabsorption cross sec- 
tion with VijCOM to those obtained with other potential 
models: the AV18 potential and the AV18 with the in- 
clusion of the Urbana IX three-body force, AV18+UIX. 
Predictions for 3 H and 4 He are taken from [32J and [lj, re- 
spectively. The addition of the short range non-locality 
introduced by the correlators in the UCOM has conse- 
quences on er 7 . In fact, with respect to the usual AV18 
potential, the Vucom leads to a reduction of the peak 
cross section of about 4% for 3 H and 10% for 4 He, while 
a general enhancement of the tail of the cross section 
is found: it amounts, for example, to 15% and 23% at 
u> = 60 MeV, for 3 H and 4 He, respectively. Similar effects 



6 



< 



< 



0.2 
0.1 



-0.1 
-0.2 
-0.3 

0.2 
0.1 



-0.1 
-0.2 
-0.3 




20 40 60 80 100 120 140 

uj [MeVl 





4 He - 


\ 




\ 









40 60 80 100 120 140 

uj [MeV] 



FIG. 5: (Color online) Quantity A potcntial (tj) (see text) as a 
function of the photon energy in case of 3 H and 4 He for the 
Vucom and AV18 potentials. 



on the photoabsorption cross section are obtained with 
the introduction of the UIX three-body force, though 
with some differences. With respect to the interaction 
model of AV18+UIX the Vucom potential leads to a 4% 
higher photodisintegration peak in case of 3 H and to a 
4% lower one for 4 He. We also observe a 0.5 and 2 MeV 
shift of the peak position towards higher energies for 3 H 
and 4 He, respectively. At low photon energies predic- 
tion of Vucom are lower than AV18+UIX: for example 
the difference is about 10% for 3 H and 25% for 4 He at 
3.5 MeV after disintegration threshold. Interestingly, the 
tail of the photoabsorption cross section obtained with 
the Vucom is very similar (less than 6 — 7% difference) 
to the AV18+UIX result for photon energy 20 < u < 120 
MeV in case of 3 H and 45 < u < 100 MeV for the alpha 
particle. 

We observe that with the Vucom a reduction of the 
missing genuine three-body force is achieved for high pho- 
ton energies. In fact, keeping as reference the AV18+UIX 
curve, the relative difference on er 7 , defined as 



y\ potential (f,,^ 



1 



—potential/. ,\ 



AV18+UI 



'1 



(11) 



gives information on that. In Fig. [5j the quantity 
AP° tontial (w) is presented in case of 3 H and 4 He for uj 
beyond the corresponding peak energies obtained with 
Vucom- In case of 3 H, AP°tentiai( w ) is reduced by at 
least a factor 2 in the energy region 25 < uj < 75 MeV 
and 90 < uj < 120, going from the AV18 to the Vucom 



Potential 


(rl) 3 H 


(r 2 ) 4 He 


(r 2 vv ) 4 He 


Vucom 


2.52(1) 


1.96(1) 


5.41(1) 


AV18+UIX 


2.51(1) 


2.05(1) 


5.67(1) 



TABLE II: Mean square point proton and proton-proton radii 
in fm 2 with Vucom and AV18+UIX ((r 2 ) = (r 2 ) for 4 He 
without isospin-mixing) . Results for triton from [50( and for 
alpha particle from [5lj. 



potential. An even stronger reduction effect is found for 
4 He for energies between 45 MeV and pion threshold. 

In the following an investigation of two photonuclear 
sum rules, the already mentioned m_i ( Eq. (10)) and the 
Thomas-Reiche-Kuhn (TRK) sum rule [l| 



m = 



2J, 



ct 7 (w) duj 
o + 1 Tr x 



(12) 



is briefly presented (see also Ref. [51|). With respect to 
the AV18+UIX case, the Vucom potential predicts a very 
similar m_i for 3 H and a lower m_i for 4 He. This fact 
is also reflected in the size of the radii of Eq. (fTU|) . In Ta- 
ble|TT]we summarize the situation: for 3 H we present only 
the proton radius, since the there is no proton-proton 
radius, while for 4 He we show the squared mass radius 
((r 2 ) = (r 2 ) if isospin mixing is neglected) and (fpp). 
One can see that Vucom leads almost to the same pro- 
ton radius as AV18+UIX within a 0.4% deviation, but 
smaller radii for the alpha particle are obtained, if com- 
pared to AV18+UIX. A consequence of that is the lower 
4 He peak cross section found with Vucom with respect 
to AV18+UIX. The TRK sum rule is of interest since 
it contains information on the exchange mechanisms in- 
duced by the non-commuting part of the potential V. For 
example, for 4 He, we find that too = 117.6 mb MeV for 
Vucom, which is about 20% lower in comparison to the 
value of 146.2 mb MeV found for AV18+UIX [p. This 
is more strongly connected to the fact that with Vucom 
the cross section goes faster to zero beyond pion thresh- 
old and indicates that the meson exchanges underlying 
the two potential models are different. 

We now compare the theoretical predictions with the 
available experimental data. Among the theoretical 
curves we add the result with the JISP potential [2l| . 
recently published in (48[. In Fig. [51 we show the sit- 
uation for the 3 He nucleus and for 4 He. In case of the 
three-body nucleus the error bars of the data from Feti- 
sov et al. [53j are unfortunately too big to allow us to 
discriminate among the different potential models. 

In case of the alpha particle the experimental situa- 
tion is more involved. Close to threshold several data 
were taken in different experiments, which unfortunately 
show fairly large discrepancies. Only the data from Arka- 
tov et al. j&ij cover a larger energy range. They present 
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FIG. 6: (Color online) Theoretical photoabsorption cross sec- 
tion of 3 He and 4 He with three different potential models, 
(thin line) Vucom, (thick line) AV18+UIX and (dashed line) 
JISP potential from [48|. in comparison with the available 
experimental data: empty squares fro m [531 ] for 3 He, empty 
triangles from [54j], full triangles from [Jfull squares from 

]|, shaded area sum of data from [55j and Hi and box from 

53 for 4 He. 



a rather high peak cross section, in fair agreement with 
recent data from 0, which favors the AV18+UIX poten- 
tial model. In the energy region 30 < oj < 45 MeV 
the results with Vucom agree better with this set of 
data than those with AV18+UIX. Finally, at higher en- 
ergies the measurements are precise enough to conclude 
that AV18+UIX and Vucom can explain, within the er- 
ror bars, the photoabsorption cross section, in contradis- 
tinction to the JISP potential. In Ref. [48| it was ar- 
gued that a possible reason of this discrepancy of the 
JISP model could be due to the probably incorrect long 
range part of the potential, whose construction is not 
constrained by the meson exchange theory. Regarding 
this issue we can only state that Vucom, which explic- 
itly contains the long-range pion exchange term of AV18, 
unchanged by the short range action of the correlations, 
exhibits a correct tail behavior in the photoabsorption 
cross section of light nuclei. Here we would also like to 
point out that the consideration of large Hilbert spaces 



(large K max and n max ) is essential to achieve a stable 
and HO-indcpcndcnt result on the whole photon energy 
range. The fact that the JISP potential is constructed on 
a small HO basis, for defined n max and HQ values, could 
limit the correct description of the tail of this sensitive 
observable. 



IV. CONCLUSIONS 

We have presented the results of an ab initio calcu- 
lation of the 3 H, 3 He and 4 He photodisintegration cross 
section with the Vucom potential. The difficulty of the 
scattering continuum problem is circumvented via the 
Lorentz Integral Transform method, where the final state 
interaction is fully taken into account. An expansion of 
the wave function on hypersphcrical harmonics, extended 
for non-local interaction, is used for the solution of the 
"Schrodinger-like" equation. The first investigation of 
the Unitary Correlation Operator Method on a contin- 
uum observable is presented in this paper. The sensi- 
tivity of the photoabsorption cross section on non-local 
terms of the interaction is investigated. With respect to 
the traditional local AV18 potential, a reduction of the 
peak and an enhancement of the tail of the cross section 
is produced by the non- locality of Vucom • The compari- 
son between the AV18+UIX potential model and Vucom 
allows to investigate, whether the omission of a genuine 
three-body force can be replaced by a non-local interac- 
tion on a continuum observable. Though binding energies 
of very light nuclei are well described by Vucom with a 
1 — 2% difference with respect to the AV18+UIX model, 
the situation is different in the photodisintegration cross 
section below pion threshold. Our analysis shows that 
the Vucom potential leads to a very similar result as the 
AV18+UIX for high photon energy, while in the region 
close to threshold larger differences are found, particu- 
larly in case of 4 He. The larger deviation in 4 He with 
respect to the 3 H case could be related to the fact that 
smaller mean square radii are obtained with the Vucom ■ 
The higher density of the alpha particle makes it certainly 
more sensitive to characteristic short range properties of 
the three-body force and non-local parts of the NN inter- 
action, magnifying the differences. The similarity of the 
high-energy cross section shows that also in the contin- 
uum genuine three-body forces can be partly simulated 
with the non-local Vucom potential. But the differences 
found at low-photon energies point out that the two po- 
tential models are not completely equivalent. Unfortu- 
nately, the lack of precise experimental data limits the 
possibility to use such calculations as a discriminant test 
of the potentials. 
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